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Abstract. We present model equations for the trapping and accumulation of 
particles in a short cylindrical channel with nanostructured inner walls when a fluid 
passes through, carrying a moderate load of impurities. The basic ingredient of the 
model is the introduction of a phenomenological "effective-charge density" of the walls, 
related to the electrical charges exposed in the nanotexture, and which is gradually 
reduced as the flow runs through the channel and the trapped impurities cover the 
internal walls. By solving the proposed equations, three regimes are predicted for the 
channel: a linear or clean- fllter regime, a logarithmic or half-saturation regime, and the 
saturation limit. It is proposed that experimentally testing these regimes may help to 
understand the enhanced trapping capability observed in many diverse nanotextured 
channel structures. 
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1 Background introduction 



Recently, [lj|8| researchers from both academia and industry have experimentally demon- 
strated that a variety of nanostructures and nanotextured media (see some examples 
later in this introduction) can efficaciously trap nanoimpurities carried by fluids when 
a flow is induced by external hydrostatic pressure. These flndings are not only sci- 
entiflcally interesting, but also promising for the socially and economically important 
application of puriflcation of drinking water, j4||9] and of other liquids ^ . When com- 
pared to conventional porous fllters, the new media have the important advantages of 
retaining impurities of sizes typically in the tens of nanometers and, at the same time, 
presenting a resistance to hydrodynamic flow orders of magnitude smaller than what 
conventional models would predict for channels of diameters as small as the particles 
being trapped. 

Roughly, we can divide the structures presenting such enhanced impurity trap- 
ping capability into two groups: a) Those formed by nanometric-diameter channels 
through which the fluid flows. [l}]4] A well known example is the nanotube arrays 
grown and experimentally tested by Srivastava and coworkers. |l] Other specially in- 
teresting example are graphene membranes, although by now they have been probed 
only through molecular dynamics simulations. [2] In any nanometric-diameter channel, 
simple size exclusion will play a major role in the retention of nanoimpurities. Still, 
these structures exhibit remarkable trapping capability also for some ions signiflcantly 
smaller than the channels' diameter. (l||2] The resistance to flow is observed to be 
well lower than what conventional models predict for these diameters, a phenomena 
often attributed to water-nanostructure interactions (see, e.g., jl]) though not yet fully 
understood at the quantitative calculation level, b) The second group corresponds to 
nanostructures embedded in larger structures, resulting in fllters composed by channels 
with micrometric diameters and inner walls coated with nanoparticles. Examples are 
conventional microfllters coated with Y2O3 [5], Zr02 [6], or AI2O3 [7|[8] nanopowders 
(further examples can be found in the reviews [3||4|[9]). These structures have been 
observed by their growers to have a surprisingly good flltration performance for nano- 
metric impurities, as small as ~ 10 nm, in spite of the relatively large diameter of the 
channels (note that in a channel of diameter 1 /xm only about 0.04% of the fluid will 
transit closer than 10 nm from the walls). l3]-[9] Their hydrodynamic resistance is quite 
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low, similar to the one of conventional micrometric filters. The filtration mechanism 
of these structures is unknown as yet. Nonetheless, the initial trapping capability was 
observed to depend on pH, [5]-[7] and thus electrostatic and polar attraction may be 
suspected to play a role. 

The purpose of this Letter is to propose a semi-macroscopic model for the nanoim- 
purity trapping capability of short cylindrical-like channels with nanostructured inner 
walls of the type composing filters of category b in the previous paragraph. The equa- 
tions do not assume any detailed chemical compound for the nanotexture, in corre- 
spondence with the ample variety of nanostructures producing the enhanced impurity 
trapping capability. The model produces agreement with the trapping performances 
quantitative reported by experimentalists in various systems. These performances cor- 




charmel nontrapped impurity 



Figure 1: Representation of a nanostructured channel filter as modelled in the present 
Letter. The radiuses tq and po correspond to the average dimensions of the bare 
channel and impurities. The effective radiuses Ve and pe vary as trapped impurities 
cover the inner wall, via their dependences on, respectively, the areal density n of 
trapped impurities and on the areal density Zg of effective charge of the inner wall. 
This Zf. reflects that exposed charges in a nanostructured surface attract the impurities 
in the fluid and also constitute binding anchors for those impurities. It is expected 
to diminish as impurities cover the surface, for which we assume the simple Ze{n) 
dependence plotted at the right (although it can be easily generalized to more complex 
functionalities, see main text). 
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respond always to the initial filtration behaviour, i.e., they were always measured in 
the clean-channel limit. In order to further explore the model, we also calculate the 
time-evolution of the trapping performance and impurity accumulation as flow contin- 
ues and the channels become dirtier. This allows us to propose future measurements 
that may test these ideas more thoroughly. We believe that some aspects of the model 
could be also useful to explain in part the trapping of the smaller ions in the nanodi- 
ameter channels of category a. But its full applicability to that case is limited by our 
use of classical dynamics for the carrying fluid, and hence we do not focus here on that 
category (also, for these nanodiameter channels, in which the number of fluid atoms 
is manageably small, molecular dynamics simulations as those in [2] could be a more 
reliable, albeit not general, approach). 

Our theoretical proposal explores a simplified (but realistic) view, in which the im- 
proved filtration capability is primarily due to the fact that the nanotexturing exposes 
electrical charges in the walls which induce both electrostatic and van der Waals at- 
traction over the impurities in the fluid. This nanostructuring also provides chemical 
anchors for the binding of those impurities once they collide with the channel walls. 
Correspondingly, our basic ingredients will be the introduction of an effective-charge 
density, z^, of the inner walls of the channels, and writing down as a function of Zg the 
probabilities for impurity collision and binding in the channel walls, ficoiude and fibmd- 
As it could be expected, Ze will vary with the areal density n of impurities trapped 
in the inner walls of the channel. We then obtain within the model the evolution of 
n with time t, as (for the fully explicit form of this equation, see Eq. ^ later in this 
Letter) : 



where $imp is the impurity flow from the incoming fluid. This equation is nonlinear. By 
solving it both numerically and through analytical approximations, we predict three 
different regimes on the n{t) evolution: a linear or clean-filter regime, a logarithmic or 
half-saturation regime, and the saturation limit (see Fig. 2a). 



dn 
dt 
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2 Obtainment of an equation for the areal density 
of trapped impurities in a channel with nanos- 
tructured walls 



2.1 Initial modelling and notations 

Our starting point, and most of our basic notations, are illustrated in Fig. 1. We 
consider a channel with nanostructured inner walls, its nominal shape being cylindrical- 
like with average radius tq and length 6x. Inside this channel a fluid flows due to 
externally-applied hydrostatic pressure, carrying a load of impurities. Some of those 
impurities will become trapped by the inner wall of the channel, then reducing its 
effective radius to a value r^. = Tq — rin, where n is the number areal density of trapped 
impurities in the inner wall and ri is a constant of the order of the average radius 
of impurities, po- All through this Letter, by "areal density" we refer to quantities 
normalized using the nominal area of the inner wall, 27rro6x (and not the cross-section 
of the channel). Also, for simplicity wc consider all impurities equal among them 
(subsequent generalization to multiple chemical species should be easy) and the length 
6x short enough as to take n and tq constant through the channel's axis coordinate 
(generalization to longer lengths, and also to some different shapes, should also be easy, 
by integration over dx). The impurity concentration in the ffuid is considered to be 
moderate enough as to not signiflcantly affect its viscosity, and as for the impurities in 
the fluid to be noninteracting with each other (specially when colliding with the inner 
wall) . 



2.2 Effective-charge density of the inner wall, Ze 



We now introduce the important concept of a phenomenological "effective charge" of 
the inner wall of the channel. We quantify this effective charge via its areal density 
and, as already commented on in the introduction, it reflects the fact that the 
nanostructured walls expose charges that induce both electrostatic and van dcr Waals 
attraction over the components of the impurities in the fluid. Indeed, will depend on 
the areal density of already trapped impurities n (which will screen out the wall) and 
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also on the chemistry specifics of the wall and impurities. Let us focus on the mutual 
interplays between n and Ze and in obtaining an equation for their evolution with time 
as flow passes through the channel. In particular, the interdependence Zf.{n) may be 
naturally expected to be continuously decreasing when n increases, to take a finite 
value zq aX n = {] (clean filter), and to saturate to zero when n reaches some critical 
value Tf'^^ at which all active centers of the wall become well covered by impurities. 
We thus postulate the simplest Zein) dependence fulfilling such conditions (plotted in 
Fig. 1): 



Ze{n) 



^0(1 --it), '^^n< 



n 



sat 



if n > rf^^ 



^0 1 



n 



n 



sat 



(2) 



where the notation ||...|| stands for min{l,...}. Obviously other sensible choices for 
Zein) are possible such as, e.g., Zein) =Zq{1 — with zi a positive coefficient 

that probably covers at a good approximation most actual possibilities depending on 
its value, as it corresponds to Ze{n) downwards curvature if < < 1, to no curvature 
for zi = 1 {i.e., Eq. (|2])), and to upwards curvature for zi > 1. For simplicity, we shall 
consider in this Letter only the case zi = 1. 



2.3 Impurity trapping probabilities as a function of Ze and n 

The role played by Ze in our model will be in fact two-fold. First, it affects how large 
is the distance within which if the impurity approximates the inner wall the latter 
attracts the former so much as to consider it as a collision. This attraction distance 
may be seen as an effective radius, pe, of the impurity (see Fig. 1), so that if the distance 
from the center of the impurity to the center of the channel is larger than — the 
impurity will actually touch the wall (dressed with already trapped impurities). For 
concreteness, we will assume that the dependence of pe on Zg can be linearized, in the 
relevant range of values, as Pe = Po + Pi^e, where pi is a constant. (In fact, an exact 
linear relation can be easily obtained on the grounds of simple energy balances using 
an electrostatic interaction and with the thermal energy playing the role of an escape 
kinetic energy, suggesting then also that pi will be roughly inversely proportional to 
temperature.) Using this pe we may now obtain the probability ficoiudc that, at a 
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given instant, a given flowing impurity will collide with the impurity-dressed wall. For 
that, we assume that the concentration of impurities is constant in all the transversal 



section of the fluid and that the fluid velocity profile is given by the Poiseuille law, 10 

2 



u[r) oc 



r 



where u is the fluid velocity and r the distance to the channel's 
axis (see Ref. 11 for an explicit discussion supporting that at least for channels of 
radius > 10 nm the flows of water-like liquids driven by hydrostatic pressure are in 
fact in the Poiseuille regime). The probability f^coiude is then given by the fraction 
of the fluid mass that passes through the outer ring given hj Vg — pe < r < r^, i.e., 
^collide = /rg'^-pg 'w(r)rdr/ Jq^ u{r)rdr. The result of those integrations is: 



n 
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(3) 



Note that, through the re(ra) and the pe{ze) dependences, ficoiude will vary as flow passes 
and the walls become dirtier. 



The second influence played by Ze in our model concerns the probability that an 
impurity gets actually bound to the inner wall of the channel once it actually is within 
a collision distance from that wall. The probability of such a process, per unit length, 
will be denoted here as fibmd- Obviously this probability will depend on the chemistry 
of impurities and active centers of the nanostructure, and also on the number of active 
centers not yet saturated by existing bindings. The latter indicates that fibind will 
grow with Ze, and in particular we may adopt the natural flrst-order approximation 
^bind = fio + ^i-2e (^0 rcflccts then a conventional binding probability in a conventional 
non-nanostructured fllter and VLq <^ ^i^^o)- Note that the probability fitrap that a given 
impurity in the fluid flow gets trapped in the walls during its transit trough the channel 
is then given by fitrap = ^coiudc^bind'^a^, expressed above as a function of n and Ze. 



2.4 Equation for dn/dt 



Let us now build, on the basis of the above relationships, an equation for the evolution 
of the areal density of trapped impurities, function of time when an impure 

fluid flows through the channel due to hydrostatic pressure. In fact, dn/dt itself is an 
almost equally important quantity, as it gives the flltration capability of the channel. 
This derivative may evidently be expressed as dn/dt = (<l>imp/27rro5x)fitrap, where 
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$imp is the flow of impurities brouglit by tlie incoming fluid (in units of s 
(27rro5a;)~^ is due to tlie areal density normalization in the definition of n) 
together the previous results in this Letter, this can be written now as: 



-1. 



dn 
dt 
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with Ze{n) given by Eq. Q. This Eq. Q may already be integrated if we further assume 
'^'imp(^) to be constant with n or, equivalently, with time. However, we believe that 
for the experiential applications of this equation the most interesting case will be the 
one of constant pressure difference P between both ends of the channel. Then $imp 
will depend on n and, in particular, as the channel gets dirty it may become clogged 
enough as to significantly decrease the flow of the fluid. For this case of constant P, 
the $imp(^) dependence can be obtained by employing again the Poiseuille law for fluid 
flows driven by hydrostatic pressure in cylindrical channels, that relates the channel's 
dimensions to its resistance to flow as $ = 7iPre'^/8ri6x, where $ = $imp/Cimp is the 
fluid flow, Cimp is the incoming number concentration of impurities, and rj the viscosity 
of the fluid. For the case of constant P and Cimp we then get: 
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orjox 
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(Note that for n = n^^°^ = tq/vi it is $imp 
fully closed by impurities at that value). 



and Te = 0, i.e., the channel becomes 



While Eqs. (|2|) to Q form now a closed set of equations that can be already solved 
to determine n{t), it will be very useful to flrst reexpress them more compactly, and 
in terms of a more recognizable set of physical quantities. In particular, and as it 
could have been expected, it turns out that the flltering behaviour can be described in 
terms of its features at two speciflc instants of its evolution: Its clean point, deflned 
by n'^^^^^ = 0, and its saturation point, deflned by n = n*^^* {i.e., Ze = zl'^^ = 0). In 
particular, Eqs. ([2]) to ([s]) can be rewritten as 
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where n' = n/rf^^ and t' = t/ (IGr/fe/Cimp-P) are the so-called reduced density of 
trapped impurities and reduced time, respectively. (In fact, already a first result of our 
equations is that the quantities 77, 5x, Cimp and P only affect the impurity trapping 
behaviour in the time scale, as they enter Eq. ([6]) only through the t'-normalization.) 
The other parameters in Eq. ([6]) are the effective radiuses of channel and impurities, 
and the binding probabilities, in the clean and saturated points: 



„clean 

e 

^clean 
re 

(~) clean 
^ 'bind 
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Qsat 
^ 'bind 



Po + PiZo, 



(7) 



To — riW 
Po, 
^0. 



sat 



These 6 variables are thus equivalent to the set {ro, ri, pq, PiZq, Qq, QiZq} and together 
with rf^^ they define the n'(t') behaviour through Eq. ([6]), as will be discussed in the 
next section. (Actually, in all rigour these 7 degrees of freedom could be further reduced 
to 6 by rescahng them to a single reference length, such as po, i.e., pl^^). 



3 Discussion of the n'{t') obtained by integrating 
the model equation 

3.1 Numerical integration 

We show in Fig. 2(a) an example of the n'{t') results obtained by integrating Eq. ^ 
using some representative values for the parameters involved (and always in the case 
of constant P and Cimp, and starting from a clean initial state n'(0) = 0). The com- 
putation is lightweight and takes less than a minute in current personal computers. In 
particular, we have chosen for the impurity sizes p|^* = Po = 13A and pf*^^^ = 20po; for 
the channel radius sizes rf'^'^^ = 500po and rl^^ = 450po; for the saturation impurity 
density n^^^ = and for the binding probabilities fi^imf — 0.5/lO^po and f^btnd — 



9 



(so that we neglect conventional filtration mechanisms and focus on the effects of nanos- 
tructuring alone). We have chosen these values because they reproduce the nominal 
sizes and the initial filtration capability experimentally reported for channels coated 
with Y2O3 nanopowders |5] and also for channels coated with Zr02 nanopowders [6]. In 
both cases, channel arrays of length ~ 7 mm have been show to retain MS2 viruses (of 
radius ~ 13 nm) with an initial logarithmic reduction value of 7-log when the channels 
are still clean. Note that the initial logarithmic reduction value in a finite-length clean 
channel can be simply estimated in our model through the accumulated non-trapping 
probability as -(L/po) \og^o[l - pon^^^^lipf rf'"'' - - where L is the length 




Figure 2: (a): Results, obtained by integrating Eq. ([6]), for the areal density of trapped 
impurities in units of the saturation value, n' = n/rf^'^ (continuous line), as a function of 
time in units of the half-saturation value, t' jt'^j^ with t'j^^g defined by n'{t\i^ = 0.5. The 
parameter values used are (see main text for details): pl^^ = Po = ISA, pf^'^^ = 20po, 
rf* = 450po, rf^"""" = 500po, = 5/po^ ^Itd = 0, ^u^l = O.S/lO^po, and the 
initial value was n'{0) = 0. We also show the results obtained from the clean-regime 
approximation, Eq. ([s]) (dashed line), from the half-saturation approximation, Eq. (10) 
(dotted line), and for the saturation limit, n = n^^^ (dot-dashed line), (b): Filtration 
capability dn' /dt' (continuous line) and fluid flow $ (dashed line), normalized to their 
values at t' = and as a function of obtained again from Eq. (j6| and the 

same parameter values as in Fig. 2a. Note that the horizontal axes of these plots are 
logarithmic, so that the half-saturation regime is well longer than the clean regime. 
The hydrodynamic flow through the channel is reduced only about 30% during its 
working lifetime. 
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of such channel. With the above choice of parameter values this estimate yields the 
measured |5||6] 7-log impurity trapping. Obviously a more stringent determination of 
the parameter values for these example filters would need further measurements not 
only of the initial filtration performance but also of its detailed evolution with time. 
Also obviously, other example channels may correspond to different parameter values: 
for instance, to reproduce the 94% retention rate of MS2 viruses reported in [T] in 
channel arrays with dimensions similar to the ones in [5|[6] but with a 10wt%-Al2O3 
nanocoating, the choice f^b'nd" = 0.8/lO^po would be more appropriate. 

As show in Fig. 2(b), we have also computed, using the same parameter values 
as for Fig. 2(a), the evolution with time of both the fluid flow $ and the derivative 
dra'/dt' that gives a measure of the filtration effectiveness at each instant. Note that 
$ experiences only a moderate reduction during the filtering working lifetime, other 
of the important features observed by industry researchers [7||8] and of paramount 
importance for applications. 



3.2 Clean, half-saturation and saturation regimes 

At least three regimes can be easily noticed in Fig. 2(a): The initial behaviour or clean 
regime, at which the growth of n'{t') is approximately linear; ii) an intermediate regime, 
henceforth called half-saturation regime, where the growth of n'{t') is approximately 
logarithmic; and ii) the saturation limit in which n approaches a value n^^^ at a slow 
pace. In fact, in that figure we also show as dashed or dotted lines the simple analytic 
approximations for n'{t') that can be easily found in the first two of these regimes 
by introducing in Eq. ^ the corresponding approximations t' = and, respectively, 
n' = 1/2. This equation becomes in those approximations as follows: 



n'it') ^ A[J 



dn' Af 



(linear or clean regime t' ^ 0), 



with 



A' 



clean 
bind 



lin sat^clcan 
It I g 



(^^clean ^ ^sat _ ^clean^2 _ (^^clean^2 j^g^ 



11 



and 
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and t\^2 defined by ^'(^'1/2) = 1/2 in Eq. (j6| (no truly explicit expression for t\^2 seems 
to be feasible). In obtaining the above Eqs. (|8| to (11), we have assumed that n'{0) = 
and that pe < at t' = or t[^2- 

Unfortunately, to our knowledge no measurements exist of the time-evolution of 
the filtering efficiency of channels with nanostructured walls with a precision valid 
for a quantitative comparison with the corresponding results of our equations. We 
thus propose (see below) that such measurements should be made to further clarify 
the mechanism behind the enhanced impurity trapping capability of the chanels with 
nanostructured inner walls. 



4 Conclusions and proposals for future work 



This Letter proposes a model for the main generic features of the channels with nanos- 
tructured inner walls with respect to trapping and accumulation of impurities carried 
by fiuids. This includes, e.g., their capability to clean the fiuid from impurities of size 
well smaller than the channels' nominal radius, with comparatively small resistance to 
fiow (much smaller than in a conventional channels of radius as small as the impurities). 
The model attributes the enhanced filtration capability to the long-range attraction 
exerted by the exposed charges in the nanostructured walls, and also to their binding 
capability once the impurities actually collide with them. Both features were quantita- 
tively accounted for by means of a phenomenological "effective-charge density" of the 
nanostructured wall. The model also predicts an specific time-evolution of the trapped 



12 



impurity concentration and of the filtering capability, which was solved numerically, 
and also analytically in the so-called clean and half-saturation approximations. 

We believe that our equations could make possible some valuable future work, of 
which two specific matters seem to us more compelling: First, it would be interesting to 
check at the quantitative level the agreement with experiments of the time-evolutions 
predicted above. For that, we propose to perform time-dependent measurements made 
in controlled flow setups. We have chosen in our Eq. ^ flow constraints which seem 
appropriate for this purpose; still, the model can be tested with different setups by 
just numerically integrating Eq. ^ under the corresponding experimental constraints. 
One of the prime purposes of this Letter is in fact to suggest such measurements. 

A second interesting future work (already in progress in our research group) is the 
design of optimal geometries for the combinations of channels forming filters. Our 
model opens this possibility because of the explicit use of the Poiseuille relations, that 
allow to calculate the resistance to flow of complex associations of those channels, in 
series and/or parallel. The effective diffusivity and tortuosity of the pathways network 
are also accounted for by these equivalent- circuit analyses. 
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